Connectivity transition in the frustrated 5=1 chain revisited 
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The phase transition in the antiferromagnetic isotropic Heisenberg 5 = 1 chain with frustrating 
next-nearest neighbor couphng a is reconsidered. We identify the order parameter of the large-a 
phase as describing two intertwined strings, each possessing a usual string order. The transition has 
a topological nature determined by the change in the string connectivity. Numerical evidence from 
the DMRG results is supported by the effective theory based on soliton states. 
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In recent years, there has been a standing interest 
in low-dimensional frustrated spin systems, motivated 
by their rich behavior which is often not completely 
understood. El A fundamental example of such a system 
in one dimension is the isotropic Heisenberg spin-S* chain 
with antiferromagnetic interactions between nearest and 
next-nearest neighbors, described by the Hamiltonian 



= ^ SjSi+i -I- SjSi+2. 



(1) 



The physics of frustrated chains with half-integer S is 
relatively well understood. For small a they are gapless, 
and above the certain a ~ ac a Kosterlitz-Thouless-type 
transition into the gapped dimcrizcd phase occurs, with 
the exponentially slow opening of the gap. 

For integer- S* chains the situation is much less clear. 
At a = they are in the Haldane phase with a finite 
gap to the elementary excitations. For S = 1, the Hal- 
dane phase has the so-called stmng ordera arising due to 
the broken Z2 x Z2 symmetryQ and determined by the 
nonlocal correlator 

n'-l 

0?(n,n') = -(^°(exp ^ ^7T S^^) S?,,) , a^{x,y,z). 

j=n+l 

(2) 

The physics of the S = 1 chain is believed to be well 
captured by the Affleck-Kennedy-Lieb-Tasaki (AKLT) 
modelo which differs from the Heisenberg model by the 
additional biquadratic exchange term i(SjSi-|_i)^ in the 
Hamiltonian, and has the exactly known ground state of 
the valence bond solid (VBS) type shown in Fig. |l|a. 

Several years ago it was observedO that the Haldane 
phase in S" = 1 chain breaks down at a > ac — 0.75 
in a first order transition, which is characterized by 
a discontinuous vanishing of the string order parame- 
ter (SOP) Oi = lim|„_„/|^oo Oi (n, ".') at the transition 
point, while the gap remains finite. On the basis of sim- 
ple energetic arguments exploiting the AKLT-type varia- 
tional states the tr|a,nsition in the S = 1 chain was heuris- 
tically interpretedcl as a decoupling of a single Haldane 



chain (cf. Fi£j ^) into two subchains (cf. Fig.|^b). Later 
it was showro that in the anisotropic chain the large-a 
"double Haldane" (DH) phase persists in a finite region 
of anisotropics. 

However, the order parameter of the DH phase was 
never identified, which made the real sense of the tran- 
sition rather unclear: complete decoupling would be 
reached only for a 00, and at finite a the SOP de.- 
fined on a subchain is not the proper order parameterE 

In the present paper we show,u sing the numerical 
results from the density matrix renormalization group 
(DMRG) supported by the effective theory based on soli- 
ton states, that the transition into the DH phase corre- 
sponds to the decoupling into two intertwined substrings, 
each having a usual string order, and identify the cor- 
responding order parameter. The transition thus has a 
topological nature determined by the change in the string 
connectivity. 

We start with recalling some basic facts about the 
AKLT model. Its exact ground state can b^i^-epresented 
in the compact matrix product (MP) form:l 



'^AKLT = tr{gig2 . 
g„ = (l/V3)(a+|- 



■9n}, 



(3) 



where are the Pauli matrices in the spherical ba- 
sis, fi — 0, ±1, and arc the spin states at the n- 
th site. The AKLT state possesses perfect string order, 
0° = 4/9. Elementary excitations of the AKLT chain 
are solitons in the string order.u A soliton with = ^ 
at the n-th site is well approximated by the MP state 

1^, 71) = tr{5i52 • ■ • (5«cr^)5«+i • ■ • ffA^ } (4) 

It is more convenient to use the equivalent set of states 
|r, n) with r = {x,y,z), defined by Eq. (jij) with re- 
placed by the Pauli matrices in the Cartesian basis ar- 
Soliton states |r, n) with different n are not orthogonal. 
However, one can introduce the equivalent set of states 



|r, n) = {3|r, n - 1) + |r, n)}/2V2 



(5) 



which have the property (r, n|r', n') = Srr'Snn' ■ The state 
can be represented in the same MP form (U) with 
{gncr^) replaced by the following matrix 

/; = v/273 1 |r)„ - (z/V6) £,,7 a, |0„. (6) 

Here |r)„ are the Cartesian spin-1 states at the site n, 
{r^jjl) € {x,y,z), and 1 is the 2x2 unit matrix. One 
can check that in the soliton state \r,n) the string or- 
der correlators C^l with r' ^ r change sign when 
n gets inside the {1,1') interval, while Ol{l,l') remains 
insensitive to the presence of the soliton. 

The ground state of the Heisenberg 5 = 1 chain dif- 
fers from that of the AKLT model by the presence of 
a finite density of soliton pairs. Indeed, the action of 
the spin operator on the AKLT state is to create 
the state + |/i,n— 1). Thus generally the action 

of the Hamiltonian produces states with soliton pairs of 
the type |r, n;r, n'), and only in the special case of 
the AKLT model the contributions of bilinear and bi- 
quadratic terms cancel each other. One may say that 
the AKLT state is a skeleton state for the Haldane chain, 
which gets dressed with the soliton pairs. 

The variational energy of the VBS state of Fig. |l|a for 
the Hamiltonian (0) is —A/3 + 2a/9 (per spin). It is easy 
to seed that there is another VBS state, being a prod- 
uct of two VBS strings defined on the a-subchains, as 
shown in Fig. whose energy — 4a/3 becomes lower 
for a > 3/4. Again, the actual large-a ground state dif- 
fers from the skeleton state of Fig. ||b by the finite density 
of soliton pairs. There are two types of pairs, with the 
solitons sitting on the same a-subchain and on different 
subchains. Presence of pairs with solitons residing on dif- 
ferent a-subchains destroys the normal SOP defined on 
a single subchain. One can check that creation of such 
pairs is equivalent to adding a finite admixture of states 
with intertwined VBS strings shown in Fig. |l|c,d. 

There is, however, another order parameter surviving 
pair creation. One can define double string order as 

m-2 

0^(n,m) = (5;j5;j+i(exp ^ z^5;)5^^„i5;^), (7) 

j=n+2 

and it is easy to check that it is not sensitive to the 
presence of any type of soliton pairs. The correspond- 
ing order parameter O2 = lim|^_„|^oo C2 (''^> is fi- 
nite in any state which is a product of exactly two VBS 
strings, arbitrarily intertwined, because one can always 
factorize (|^) into a product of two normal SOPs defined 
along those strings. On the other hand, the correlator 
decays exponentially in the AKLT state with only 
one VBS string, because in this case (0) factorizes into 
(5°5^)C'i(m— l,n+ 1). One may say that O2 measures 
the connectivity of the state, telling us how many VBS 
strings are there. It is natural to assume that O2 should 
be the order parameter characterizing the DH phase. 

The above assumption is readily confirmed by the 
DMRG results. In our DMRG calculations, we have con- 
sidered chains with open boundary conditions and up to 
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FIG. 1: Schematic picture of the ground state of the frus- 
trated 5=1 chain (hi): (a) VBS state with one string corre- 
sponding to the Haldane phase; (b)-(d) VBS states with two 
strings contained in the DH (large-a) phase. 

L = 300 spins, keeping M — 400 states. We have calcu- 
lated the excitation gap and both the conventional SOP 
Oi and the new SOP O2 (see Fig. |). While the results 
were converged in the number of states kept, we carried 
out finite size extrapolation for frustrations very close to 
the transition point. We observe that both order param- 
eters exclude each other in the sense that if either of them 
decays to zero, the thermodynamic limit of the other is 
non-zero and vice versa. The value of O2 is strictly zero 
below a < 0.77, and for a = 0.775 already at 0.061, i.e. 
44% of its asymptotic value of 0.1401 for a 00, where 
one has two completely decoupled unfrustrated chains, 
and O2 is simply the square of the string order for the 
unfrustrated S = 1 Heisenberg chain, 0.3743. In fact, for 
the frustration a — 1.3, it reaches 90% of the asymptotic 
value. At both sides of the transition, both order pa- 
rameters exhibit clearly finite values, signalling the first 
order transition. The finding of a finite gap for all values 
of a is consistent with this picture. 

In order to describe the physical picture sketched 
above, we have constructed the effective theory based 
on the soliton states (||). For a < | we regard the AKLT 
state as the vacuum state, and introduce bosonic oper- 
ators ijj^ creating states \r,n), r = {x,y,z). Similarly, 
for a > I we treat the two-subchain VBS state shown 
in Fig. |l|b as the vacuum; this state can be conveniently 
represented in the following MP form: 

^AKLT2 ^ iy{GlG2 ■ ■ ■ Gn-iGn} , (8) 

where Gn — {gn 'S) 1) and Gn — {'^ ^ Qn) are 4x4 ma- 
trices, and the total number of spins N is assumed to be 
even. The operators tjj^ are interpreted as creating soli- 
ton states \ryn) on the subchains, defined by replacing the 
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FIG. 2: (a) the string order d and double string order O2 as 
functions of the frustration q; (b) the same for the excitation 
gap A. Symbols correspond to the DMRG data, and solid 
lines are theoretical predictions from the effective model. 



The Hamiltonian does not take into account any in- 
teraction between the solitons. The most important con- 
tribution to the interaction comes from the constraint 
that at most one particle can be present at a given 
site, which is effectively equivalent to the infinite on- 
site repulsion U . We treat the effect of the constraint 
using Brueckner's appmximation along the lines pro- 
posed by Kotov et al.EJ. In this approach, one ne- 
glects the contribution of anomalous Green's functions 
and obtains in the limit U 00 the vertex function 
Trr'.ss' — ^{k,uj){SrsSr's' + Srs'^r's), whcrc fc and kuj are 
respectively the total momentum and energy of the in- 
coming particles, with 



1 



T{k,w) N 



-E 



ZqZk-qUqU^_^ 



The corresponding normal self-energy Y,{k,Lu) is 
= (4/7V) ^ Z,z;2r(fc + 9, ^ ~ f^,). 



The normal Green function has the form 



(12) 



(13) 



Gik,iv) = 



where T,± = i{E(fc, w) ± S(— fc, — w)}, and its quasipar- 
ticle part is given by 



G(fc,w) 



Zkul 



Zkvl 



l£ 



Qk - ie 



(15) 



Ti-th matrix Qn or C/„ in (^ with the matrix JF^ = (/^(g)!) 
or = [t® /^), respectively. 

Calculating all matrix elements corresponding to the 
hopping of solitons and creation of pairs, and passing 
to the momentum representation, one obtains on the 
quadratic level the effective Hamiltonian of the form 



Heff = E ^fe4-,.*fe,. + \MAA,r + ll-C-)- 



(9) 



kr 



Here for the Haldane phase (a < |) the amplitudes Ak^ 
Bk are given by the expressions 

Bk = 2(112a/27- l)cosfc 

- (10a/9)cos2fc- (13a- 3)C^^'^\ 

Ak = Bk + (4/27) X (10) 
X { 15 — 13q; + (9 — 16a) cos k — 3a cos 2fc} , 

and for the DH phase (a > |) one has 



Bk 



(2/9) (2 cos fc - a cos 2fc) - C, 



(3,2) 



ac: 



(4,2) 



Ak = Bfe + (4a/9)(5 + 3cos2fc), 
where the following shorthand notation is used: 



(11) 



^I'j) 8 3cos(/'fc) + cos((/-r)fc) 
~ 9 5-l-3cos(;fc) 



which defines the renormalization factors Zk, the Bogoli- 
ubov coefficients Uk, Vk and the spectrum flk as follows: 



rik — s_ 



Ek 



Ek = {{Ak + ^+r-Biy/', 



4 = -{l + {Ak + J:+)/Ek}, 



2 2 1 

Vk = Uk - 1, 



{Ak + S+) 
du) Ek dijj 



(16) 



where Yi± and their derivatives are understood to be 
taken sXoj — Vtk- The system of equations ([T^), ([l^), ( [l^ ) 
has to be solved self-consistently with respect to Z and 
S. This approach is valid as long as the soliton density 
P = jf ZqVq remains small, ensuring that the contri- 
bution of anomalous Green's functions is irrelevant0 

It should be remarked that our Eqs. ([l6Ldiffer from the 
corresponding expressions of Kotov et al^t3 which can be 
obtained from (nm assuming that S(fc, uj) is almost linear 
in UJ in the frequency interval {—flk,^k)] however, this 
latter assumption fails for the present model. 

As far as Eqs. (|l|), (|l|), (|l|) are solved, one can cal- 
culate the reduction of the string order caused by the 
presence of pairs. For the usual SOP in the Haldane 
phase (a < |) one finds 



Ol = {A/9){l-2jppRf 



(17) 
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Here 7 = | is the correcting factor which takes into ac- 
count that (zz) pairs do not affect Of, and the total 

density of sohton pairs Pp = ^ S«„'r(4r*i'r)^ a-iid the 
mean size of the pair R are given by 

q 

,„2 ' W^^ —}_^ZgUqVge'l .(18) 

A similar calculation for O2 in the DH phase yields 

O^, = (4/9)^(1 - jp'pR'}\l ~ 27P*i?*}', (19) 

where p^ and p* are densities of the "longitudinal" (sit- 
ting on a single subchain) and "transversal" (intersub- 
chain) pairs, respectively, 

Pp' = (Pp±Pp)/2, 
q 

and i?', Rt are the corresponding average pair sizes (note 
that all distances are defined in terms of the initial chain) : 

Em>i(2n^)«^L „t Em>o(2"^)wL+l 
R ^ —2 ' ^ = Zp • (21) 

In Fig. H we show the results of the DMRG calcula- 
tions together with the theoretical curves obtained on 
the basis of the effective theory. Though the theoretical 



calculations are not quantitatively satisfactory, they nev- 
ertheless capture the essential behavior of the system. In 
the vicinity of the transition the theoretical results are 
rather far from the numerical data; one reason is that at 
a = I the completely dimerized state has the same vari- 
ational energy as the VBS states in Fig. |l|(a,b), which is 
not taken into account in the theory. 

It should be mentioned that ajdefinition similar to (|^) 
was recently used by Todo et al.Ell for the spin-1 ladder, 
with (n, 1) and (m— 1, m) placed on two rungs. They 
have shown that such an order is present in the S — 1 
ladder for any ratio of the rung and leg exchange. One 
may thus expect that in the frustrated 5 = 1 chain with 
alternating nearest-neighbor interaction the DH phase is 
smoothly connected to the dimerized phase. However, 
this problem is not clear since in a ladder there are three 
ways to define the double string order, depending on 
whether the initial and final pairs of points are placed 
on the rungs or diagonals, and only one of the definitions 
gives a finite value in the completely dimerized state. 

In summary, we have identified the nature of the first 
order transition in the frustrated 5 = 1 chain as a change 
in the connectivity of underlying VBS states, and estab- 
lished the proper order parameter for the large-a phase. 
We have also developed the effective description based on 
soliton states, which qualitatively describes the physics 
of both phases. 
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